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Abstract 

A nonlinear Schrodinger equation, that had been obtained within 
the context of the maximum uncertainty principle, has the form of 
a difference-differential equation and exhibits some interesting prop- 
erties. Here we discuss that equation in the regime where the non- 
linearity length scale is large compared to the deBroglie wavelength; 
just as in the perturbative regime, the equation again displays some 
universality. We also briefly discuss stationary solutions to a naturally 
induced discretisation of that equation. 



1 Introduction 



Physicists study nonUnear Schrodinger equations for two purposes: (i) as 
effective equations in such fields as condensed matter pQ and (ii) to probe a 
potential deviation from exact linearity in quantum mechanics |2] ■ 

The manner in which the nonlinear Schrofinger equations are arrived at 
in the literature is varied. Here we will focus on the nonlinear equation that 
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was obtained using information theoretic arguments [3]. Such an approach, 
often termed the "maximum entropy method" or more generally "maximum 
uncertainty method" , is often used in statistical mechanics to deduce the form 
of the probability distribution. The philosophy is that one should choose 
the probability distributions with minimum bias while satisfying relevant 
constraints [1]. 

The maximum uncertainty approach had been used to understand the 
usual linear Schrodinger equation [5]. Of course in such an approach one 
must motivate the information, or inverse uncertainty, measure that is used 
to quantify the entropy (uncertainty) associated with the probability distri- 
bution. Just as the Gibbs-Shannon measure is the simplest one that satisfies 
certain axioms in statistical mechanics, the Fisher measure emerges as the 
simplest one that satifies axioms appropriate for classical ensemble dynam- 
ics, and hence for deducing non-relativistic quantum mechanics as a single 
parameter extension of classical ensemble dynamics [6]. 

One could use the approach of [6J to obtain generalisations of the Fisher 
measure, involving higher derivatives, that satisfy the same axioms but which 
involve more parameters. However the natural constraint that the equation 
be invariant under a scaling of the wavefunction then leads to nonpolynomial 
terms ~ dp/p, with p{x) the probability density; and since the nonlinearity 
is expected to be phenmenologically small, the nonlinear corrections would 
lead generally to ill-defined singular terms in perturbation theory. 

A different approach to generalising the Fisher measure is to choose the 
simplest nonsingular measure that contains the Fisher measure in some limit, 
satisfies the various axioms and which also has some physical motivation. 
That is how the regularised Kullback-Lieber measure came to be used and 
the associated nonlinear Schrodinger equation derived [3]. 

For simplicity we consider here only a single quantum particle in one 
space dimension, see ^ for the general case. The nonlinear equation is 
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and p{x) = 4'*{x)i{j{x); also 
and 

Q 

Let 

p±{x) = p{x±riL), (6) 

where the dimensionless parameter rj takes values < r/ < 1, then 

P _|_ ^ _ (1 - r])p r]P- 

{1 — ri)p + rip^ {1 — ri)p + rip^ {1 — ri)p^ + rjp 

(7) 

In addition to regulating potential singularities, t] may also be viewed as 
parametrising a family of nonlinear theories, with the usual linear quantum 
mechanics recovered at rj = 0. 

The nonlinear Schrodinger equation ([T]) resembles a differential-difference 
equation as the evolution of the wavefunction depends, any fixed time, not 
just on knowledge at the point x but also at neighbouring points a finite 
distance away, x ± riL. The nonlinearity is also non-polynomial. Neverthe- 
less, the equation shares a number of important properties with the linear 
Schrodinger equation, such as the conservation of probability and existence 
of the usual plane wave solutions. 
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2 Perturbative Regime 

It was found [7j that for smooth external potentials V{x), treating the non- 
linearity as a perturbation results in energy shifts which depend significantly 
on whether or not the unperturbed states have nodes. Those that do, have 
their energies shifted by a larger amount 

^ - V) (1 - 4r7) E + 0{L/a)\ (8) 

where a is the characteristic length scale of the linear theory, and L/a is 
the perturbation parameter. The coefficients C„p depend on the slope of 
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the unperturbed wavefunctions near the nodes, and hence on the external 
potential. Notice that the dependence of the energy shifts on r] is universal, 
resulting in positive energy shifts for small rj and negative shifts for larger 
values. 

One may also try to fix the value of [8]. If one assumes to be a 
parameter that "fiows" in the space of theories, then one may see whether 
there is some value that minimises the energy of the system. Indeed, the 
energy shifts mentioned above reach a unique global minimum at 

r]m ~ 0.80 . (9) 

At this value of t] the leading energy shifts are negative so that the nonlin- 
earity reduces the energy of the original linear system. 



3 Nonperturbative regime 

For L large, p{x ± L) ~ p{L) ~ 0, and so the time-independent nonlinear 
Schrodinger equation for an eigenstate reduces to 

E(f){x) = - ^ log(l - r])(j){x) , (10) 

where V{x) is a given external potential. This equation suggests that 
can be nonzero only at those points where the potential vanishes. For a 
potential with a single minimum, this means that the wavefunctions collapse 
to a delta function and the energy is given in this limit by 

E = -^l0g(l-77) (11) 

As a function of rj this energy attains a minimum at the universal, that is 
V{x) independent, value 

r/~0.9. (12) 

Notice also that because of (jlj), the energy goes to zero as L — > cxd. 

A variational calculation gives a quick, though approximate, interpolation 
between the small and large L regimes. We did this for the SHO potential 
using trial states which are SHO stationary states but with the usual de- 
Broglie wavelength a = ^ ^ replaced by a variational parameter b |9] . The 
dimensionless quantity c = b/a was varied to minimise the action, which 
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depends on the free parameters t] and e = L/a. The energy functional was 
then evaluated at the minimum for each [i], e). 

For small e the results are in agreement with the perturbative results as 
expected. For larger e we found that the wavefunction shrinks and expands 
suddenly at certain [r] dependent) values of e before eventually shrinking 
(c ^ 0) to a highly localised state as suggested by flTUl) . The oscillatory 
behaviour of the wavefunction might be an artifact of the variational calcu- 
lation. Figures (1-3) summarise some of the variational results. (The depen- 
dence of the energy on r] is qualitatively similar to that in the perturbative 
regime) . 



4 Discretised States 

Since the nonlinear equation already has the form of a difference-differential 
equation, it is most natural to consider the fully discretised equation, with 
p{x + L) = Pn+i etc., as another approach to study nonperturbative aspects 
of the equation [TU] . We remark that if L is related to gravity, as speculated 
in [3], then this discretisation of space suggested by the structure of the 
equation is consistent with discretisations implied by other approaches to 
quantum gravity. 

Consider first the the time-independent nonlinear equation for rj = 1. We 
have 



E = £ 



m +1 

Pn+l Pn 



(13) 



If we seek solutions on the full infinite line for which p„ is bounded for 
all ra, then by using the iterative nature of f|T3|) one easily shows that that 
this is only possible for E > and even then there is a lower bound on n (to 
keep the probability positive). That is, the discretised space is necessarily 
truncated. Similar results hold for the general regularised case, r/ 7^ 1, and 
also for the q— deformed equation in Ref.[3]. 



5 Summary and Outlook 

The fact that the nonlinear equation ([T]) was motivated through a philosoph- 
ically appealing information theoretic approach suggests that the equation 
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might eventually be of physical relevance, for example as an effective equa- 
tion in one or both of the domains mentioned in the first paragraph of the 
introduction. As such, it is of interest to investigate the properties of that 
equation. 

We found that for large e = L/a the nonlinearity acts as a strongly at- 
tractive force, causing a localisation of the wavef unctions. The resultant 
energy of all eigenstates is minimised at a universal value of ~ 0.9, which 
is larger than the value rj ~ 0.8 found (for states with nodes) in the per- 
turbative regime of e. A natural question is whether for any L, there is a 
V(x)-independent rj value for which the energy of eigenstates is minimised. 

There are several other properties of the nonlinear equation that have yet 
to be studied. Of particular interest are its time-dependent states: dynamics 
of wavepackets and possible solitary waves. We invite interested readers to 
the adventure. 

Acknowledgement R.P. thanks the organisers for their hospitalty and 
for staging a stimulating event. 

References 

[1] C. Sulem and P. Sulem The nonlinear Schrodinger equation:self-focusing 
and wave collapse (Springer, 1999). 



[2] For some reviews, see G. Svetlichny, arXiv:quant-ph/0410036 R. Car- 



roll, arXiv:quant-ph/0401082 and references therein. 



[3] R. Parwani, Ann. Phys. 315, 419 (2005). 

[4] E.T. Jaynes, Phys. Rev. 106, 620 (1957); 108, 171 (1957); 

Probability Theory, The Logic of Science (Cambridge University Press, 
2004). 

[5] M. Reginatto, Phys. Rev. A58, 1775 (1998); Erratum ibid. A60 1730 
(1999); 

B. R. Frieden, J. Mod. Opt 35,1297 (1988); Am. J. Phys. 57 (1989) 1004. 



[6] R. Parwani, J. Phys. A:Math. Gen. 38, 6231 (2005); 
R. Parwani, Int. J. Theor. Phys. 45, 1901 (2006). 



6 



[7] R. Parwani and G. Tabia, J. Phys. A: Math. Theor. 40 5621-5635 (2007). 

[8] R. Parwani, Theoretical and Mathematical Physics, 152(1): 10121016 
(2007). 

[9] R. Hasson and D. Richards, arXiv: quant- ph / 001 2068| v^l . 
[10] R. Parwani and H. S. Tan, Phys.Lett.A363:197-201 (2007). 



7 



Figure Captions 



• Figure 1 : The variationally computed energy versus e. Plots are similar 
for the ground state and the example of an excited state, n = 5, studied 
here. Generally, for fixed e, the energy decreases as r) increases. 

• Figure 2 : The normalised width c of the variational ground state versus 
e. Looking at the intercept on the x-axis from the right, the curves cor- 
respond to the following values of rj : 0.1, 0.2, 0.999999, 0.5, 0.999, 0.9. 

• Figure 3 : The normalised width c of the variational excited (n = 5) 
state versus e. 
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Figure 1: . 
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Figure 2: . 
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Figure 3: . 
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